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Abstract 

This paper develops the geometry and analysis of the averaged Euler equa- 
tions for ideal incompressible flow in domains in Euclidean space and on Rie- 
mannian manifolds, possibly with boundary. The averaged Euler equations 
involve a parameter a; one interpretation is that they are obtained by en- 
semble averaging the Euler equations in Lagrangian representation over rapid 
fluctuations whose amplitudes are of order a. The particle flows associated 
with these equations are shown to be geodesies on a suitable group of volume 
preserving diffeomorphisms, just as with the Euler equations themselves (ac- 
cording to Arnold's theorem), but with respect to a right invariant H 1 metric 
instead of the L 2 metric. The equations are also equivalent to those for a cer- 
tain second grade fluid. Additional properties of the Euler equations, such as 
smoothness of the geodesic spray (the Ebin-Marsden theorem) are also shown 
to hold. Using this nonlinear analysis framework, the limit of zero viscosity for 
the corresponding viscous equations is shown to be a regular limit, even in the 
presence of boundaries. 
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1 Introduction 

More than twenty-five years ago, using a setting introduced by Arnold [1966] , Ebin 
and Marsden [1970] proved that on manifolds with no boundary (such as spatially 
periodic flow in Euclidean three space M 3 ), the solutions of the Navier-Stokes equa- 
tion converge to those of the Euler equations as the viscosity tends to zero. Marsden, 
Ebin, and Fischer [1972] conjectured that although in a region with boundary, solu- 
tions of the Navier-Stokes equations would not in general converge to the solutions 
of the Euler equations, a certain averaged quantity of the flow may converge. Kato 
[1984] showed that the problem in the context of weak solutions of the Navier-Stokes 
equations has fundamental difficulties. 

Recently, Barenblatt and Chorin [1998a, b] also speculated that certain average 
properties of the flow possess well-defined limits as the viscosity tends to zero. One 
of the main purposes of this paper is to develop extensions of the tools of Arnold, 
Ebin and Marsden and to use them to prove that, in a sense that will be made 
precise, viscous flow, with an appropriate ensemble averaging over rapid fluctua- 
tions or spatial averaging over small scales, does indeed converge to solutions of the 
corresponding inviscid limit equations, even in the presence of boundaries. 



Some History. The geometric approach to fluid mechanics has a long and com- 
plex history, going back at least to the basic work of Poincare (see the bibliograph- 
ical references, but especially [1901b, 1910]) and the thesis of Ehrenfest (see Klein 
[1970]). In more recent times, it was the work of Arnold [1966] that was a critical 
contribution. Arnold showed, amongst other things, that if u(x, t) is a time depen- 
dent divergence free vector field on a compact Riemannian n-manifold M, possibly 
with boundary, if u is parallel to the boundary and if n(x, t) is its volume preserving 
flow, then u satisfies the Euler equations 

du _ 

— + V u u = -grad p 
at 

if and only if the curve t i— > i](-,t) is an I? geodesic in T>AM), the group of C°° 
volume preserving diffeomorphisms of M. Of course, the Euler equations in Eu- 
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clidean space in coordinates are given as follows (using the summation convention 
for repeated indices): 



du l j du l dp 
dt dxi dx % ' 

and on a Riemannian manifold, the Euler equations take the following form in 
coordinates: 

du l du l ,■ ju dp 

where gij is the Riemannian structure, g 13 is the inverse metric tensor and T l - k are 
the associated Christoffel symbols. 

Ebin and Marsden [1970] proved a remarkable result: the geodesic spray of 
the L 2 right invariant metric on V S JM), the group of volume preserving Sobolev 
H s diffeomorphisms (s > (n/2) + 1) is C°° (recall that n is the dimension of the 
underlying manifold M). They derived a number of interesting consequences from 
this result, including theorems on the convergence of solutions of the Navier-Stokes 
equations to solutions of the Euler equations as the viscosity goes to zero when M 
has no boundary. Since that time, several papers (which we shall not review here) 
reproved and in some cases, rediscovered, these results using more traditional PDE 
methods in Eulerian representation. However, the depth and beauty of the direct 
approach in Lagrangian representation remains compelling. 

There have been, of course, many other developments since that time, such 
as developments in the theory of hydrodynamic stability, bifurcations, and other 
aspects of the geometric approach. We refer to Arnold and Khesin [1998], Holm, 
Marsden, Ratiu and Weinstein [1985], Marsden [1992] and Marsden and Ratiu [1999] 
and references therein for stability studies and to Misiolek [1993, 1996] and Shkoller 
[1998] and references therein for the Riemannian geometry developments of the 
group of diffeomorphisms. 

^-Geodesies. The (k = 0) Camassa-Holm equation in one spatial dimension is 
given by 

ut - u xxt = -3uu x + 2u x u xx + uu xxx , (1.1) 



or equivalently, 



dv n 

— + uv x + 2vu x = 

dt 



where v = u — u xx . This equation has attracted much attention; amongst many 
other features, it is a completely integrable bi-Hamiltonian system with non-smooth 
solitons that have very interesting geometry. See, for example, Camassa and Holm 
[1993], Camassa, Holm and Hyman [1994], and Alber, Camassa, Holm and Marsden 
[1994]. This one dimensional equation was given an Arnold-type interpretation in 
terms of H l geodesies on an infinite dimensional group by Misiolek [1998] (who also 
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considered the case k / and used the Bott-Virasoro group) and Kouranbaeva 
[1999] (who explored the geometry of the k = equations using the diffeomorphism 
group) and Shkoller [1998] for the well posedness in H 3 / 2+e (S 1 ) (See also Theorem 
|3J3] of the present paper). 



The Averaged Euler Equations. Using recent developments in Lagrangian re- 
duction and the Euler-Poincare equations, Holm, Marsden and Ratiu [1998a, b,c], 
developed the "Arnold view" of fluid mechanics further by applying it to many 
other types of fluid equations, such as those in geophysics. One of the equations 
they studied is quite remarkable and is the subject of this work. 

The equations of concern to us may be described in two mathematically equiva- 
lent ways. This equivalence is the .ff ^analogue of the theorem of Arnold mentioned 
above: 

1. First of all, of course, there are the explicit averaged Euler equations (these 
are written in Euclidean coordinates on Euclidean n-space; we shall write the 
equations on manifolds in the main text): 



dv a dv 2 
— — h U J — — r — a 
at oxi 



dx' 1 



Au, = -^:, (1.2) 



where a is a positive constant, v = u—a 2 Au and A denotes the componentwise 
Laplacian, and there is a summation over repeated indices (and in Euclidean 
coordinates, there is no difference between indices up or down). 

As with the usual Euler equations, the function p is determined from the con- 
dition of incompressibility: div u = 0; the pressure satisfies a Poisson equation 
that is determined by taking the divergence of the equation. 

The no slip boundary conditions u = on the boundary will be assumed for 



equation (1.2). A second possible choice will be briefly discussed later. 



2. The second description is to say that the flow rjt(-) '■= vi'ify °f the time 
dependent vector field u is a geodesic in a subgroup of T> S (M) with respect to 
the right invariant H 1 metric. 

There are two interesting and distinct ways of deriving the equations from the point 
of view of mechanics: 

1. The first derivation is by means of averaging the Euler equations over rapid 
fluctuations whose amplitude is of order a (see Holm, Marsden and Ratiu 
[1998b] for the case of Euclidean space). This method is explored in Marsden, 
Ratiu and Shkoller [1999b]. 

2. The second derivation uses the notion of second grade fluids, as will be dis- 
cussed in the last section of the paper. 
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The Main Results of the Paper. In this paper we show that the geodesic spray 
of the right invariant H 1 metric on compact manifolds with boundary is smooth 
and has a unique flow which corresponds to the solution of the averaged Euler equa- 
tions (|l,2|), the correspondence given by right translation. This has consequences 
that are similar to those for the usual Euler equations, namely, local existence 
and uniqueness, smooth dependence on initial condition and time, convergence of 
Chorin-Marsden type product formulas, etc. 

Using this set up, we establish one of the main corollaries of our results; namely, 



for the viscous analog of the averaged Euler equations (|1.2|) , 

dv? 



dv 1 . ,■ < dv 1 9 



dx" 1 



which we call the averaged Navier-Stokes equations, or the Navier-Stokes-a equa- 
tions, we prove that 

The solutions for the corresponding viscous problem converge to those 
for the ideal problem, as the viscosity goes to zero, even in the presence 
of boundaries on uniform time intervals [0, T], for T > 0, independent 
of the viscosity. The size of the interval [0, T] is governed by the time of 
existence for the averaged Euler equations with the initial data fixed. 

The inclusion of boundaries is a major difference from the situation with the usual 
Navier-Stokes equations and the Euler equations, for which such convergence is 
believed to not hold near boundaries because of the generation of vorticity at the 
boundary. (See, for example, Marsden, Ebin and Fischer [1972] and Chorin, Hughes, 
Marsden, and McCracken [1978] for discussions). 



Product Formulas. Once one has the above geometric setting and the smooth- 
ness of the spray is established, one can "read off" a number of interesting conse- 
quences. One of these is a product formula stated in the following equation. This is 
one of several possible product formulas, some of which are useful in computational 
settings. The product formula we have in mind is the mean motion version of a 
formula for the Euler equations having its origins in Chorin [1969] and stated, along 
with a very simple proof (based on the smoothness of the Lagrangian representation 
of the equations) in Marsden, Ebin and Fischer [1972]. This formula is 

Et= lim {P a oG t/n ) n , 

where Et is the flow on the space of divergence free velocity fields u of the averaged 
Euler equations (with, say, zero boundary conditions), P a is the .ff 2 -orthogonal 
projection onto the divergence free vector fields zero on the boundary (defined be- 
low), and where Gt is the unconstrained H 1 spray — that is, the problem with the 
incompressibility condition dropped. 

If there is no boundary, then P a is the same as P, the usual Hodge L 2 orthogonal 
projection. For the case of zero slip boundary conditions, we have 

P a = (1 — a 2 A)~ 1 P(l — a 2 A) 
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where the domain of (1 - a 2 A) is H 2 (TM) n Hq(TM) (the Sobolev space of H 2 
vector fields vanishing on the boundary). Notice that P a is idempotent, so that it 
is a projection in H 2 since 

((1 - a 2 A)P a u, (1 - q 2 A)n)o > , 

where (•, -)o is the L 2 inner product. Note, however, that the H 1 inner product of 
(1 — a 2 A) _1 grad p with divergence free vectors U in H s n Hq vanishes in the sense 
that 

((1 - a 2 A)~ 1 grad p, U) x := ((1 - a 2 A)(l - a 2 A) _1 grad p, U) = (grad p, U) = 0. 

In the case of the viscous version of the equations, it is very interesting that now 
one does not require the vorticity creation operator to correct the boundary term 
(see also Marsden [1973, 1974], Chorin, Hughes, Marsden and McCracken [1978] 
and Benfatto and Pulvirenti [1986]). The form of this product formula is 

F t = lim (S t/n o E t/n ) n 
= lim o (5 t/n oP%G i/n ) n 

where St is the Stokes-a flow. 

Smoothness of the Spray. A crucial central result for this paper is the smooth- 
ness of the spray of the averaged Euler equations. Concretely, when we refer to the 
spray, we mean the infinite dimensional vector field governing the dynamics written 
in Lagrangian representation. It is this vector field, which is shown to be smooth in 
the standard sense of smooth vector fields on infinite dimensional manifolds. 

Smoothness of the spray has the direct consequence that the averaged Euler 
equations are well-posed using the dynamical variable u, for short time in the H s 
topology for s > n/2 + 1. Note that the vector field v = (1 — a 2 A)u is only in 
class H s ~ 2 and need not be C°. Smoothness also plays a fundamental role for the 
derivation of the consequences already mentioned, such as the limit of zero viscosity 
and the product formulas. 

Some additional results which are a consequence of the smoothness of the spray 
are also worth stating. For example, if 7/ is a volume preserving diffeomorphism close 
to the identity, then there is a unique initial velocity field whose averaged Euler flow 
reaches r\ in time one (and stays in a neighborhood of the identity diffeomorphism). 
In addition, even if the initial condition is only H s , the flow is infinitely differentiable 
in time. This latter fact is surprising even for the usual Euler equations.^] 

We conclude by outlining the general strategy for the proof of smoothness of this 
spray as well as stating a number of related results that follow from our method of 
proof. 

L As was pointed out by T. Kato (in an unpublished manuscript), this fact for the Euler equations 
was rediscovered by several authors who were evidently unaware that it is an immediate consequence 
of Ebin and Marsden [1970]. 



6 



1. The spray of the Camassa-Holm equation on the circle is smooth on the tangent 
bundle of the H s diffeomorphism group of the circle, for s > 3/2. This follows 
from the fact that the spray has no loss of derivatives (see Remark 3.5 of 
Shkoller [1998] and Kouranbaeva [1999]) and the method of proof in the present 
paper. We also note that the noncompact case (replace the circle with the real 
line) should follow using weighted Sobolev spaces, as in Cantor [1975]. 

2. More generally, the right invariant H l metric on the full H s diffeomorphism 
group for s > n/2 + 2 of a compact Riemannian n-manifold possibly with 
boundary and with pointwise fixed boundary conditions, has a smooth spray. 
Using the Moser [1966] estimates on products, we conjecture that the condition 
on s can be weakened to s > n/2 + 1. 

3. The right translated projection P a as a map from the tangent bundle of the 
full diffeomorphism group (restricted to the volume preserving ones) to the 
tangent bundle of the volume preserving diffeomorphism group is smooth. 

4. The sprays on the group of volume preserving diffeomorphisms and that on 
the full group of diffeomorphism are related by the derivative of the right 
translated projection P a . 

5. Combining the last two facts, one arrives at smoothness of the spray of the 
H 1 metric on the volume preserving diffeomorphism group in the H s topology, 
and a careful examination of the proof shows that s>n/2 + lis sufficient. 

2 Diffeomorphism Groups 

In this section we set up the relevant groups of diffeomorphisms that we shall need 
to study the averaged Euler equations in Lagrangian representation. 

Sobolev Spaces of Mappings. Let (M, (■, •)) be a compact oriented C°° n- 
dimensional Riemannian manifold possibly with boundary, and let (Q, (•, •) < ^) be 
a p-dimensional compact Riemannian manifold without boundary. By Sobolev's 
embedding theorem, when s > n/2 + k, the set of Sobolev mappings H S (M,Q) is 
a subset of C k (M,Q) with continuous inclusion, and so for s > n/2, an ilP-map 
of M into Q is pointwise well-defined. Mappings in the space H S (M,Q) are those 
whose first s distributional derivatives are square integrable in any system of charts 
covering the two manifolds. 

For s > n/2, it is known that the space H S (M, Q) is a C°° differentiable Hilbert 
manifold (see Palais [1968], Ebin and Marsden [1970] and references therein). Let 
exp: TQ — ► Q be the exponential mapping associated with (-,-)^. Then for each 
(j) G H S (M,Q), the map u cxp : T^H^M^Q) -> H S (M,Q) is used to provide a 
differentiable structure which is independent of the chosen metric, where u> cxp (v) = 
exp ov. 
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Diffeomorphism Groups. The set of H s mappings from M to itself is not a 
smooth manifold; however, if we embed M in its double M, then the set H S (M, M) is 
a C°° Hilbert manifold, and for s > n/2+1, we may form the set V S (M) consisting of 
H s maps r\ mapping M to M with H s inverses. This space is a smooth manifold. It is 
a well-known fact that the diffeomorphism group V s (M) is a C°° topological group 
for which the left translation operators are continuous and the right translation 
operators are smooth (Ebin and Marsden [1970] and references therein). One also 
knows that 77 : M — > M has an extension to an element of (the connected component 
of the identity of) V s (M) if and only if 77 lies in (the connected component of the 
identity of) V s (AT). 

Let fj.be the Riemannian volume form on M, and denote by 

V^{M) := { V eV s (M)\ V *(fi) = fi} 

the subgroup of V s (M) consisting of all volume preserving diffeomorphisms of class 
H s . For each 77 6 T> S (M), we may use the L 2 Hodge decomposition to define the 
projection P v : T V V S (M) -» T V V S ^(M) given by 

P v (X) = (P e (Xor l - 1 ))or l , 

where X G T V T)^(M), and P e is the L 2 orthogonal projection onto the divergence- 
free vector fields on M. Recall that this projection is given by 

P e (v) = v — gradp(f) — grad&(t>), 

where p is the solution of the boundary value problem 

Ap(v) = div v in M 
p(v) = on dM, 

and b solves 

Ab(v) =0 in M 

(grad b(v), n) = (v — gradp,n) on dM, 

where n is the orientation preserving normal vector field on dM. The function p 
is the pressure associated with v, while the function b is a smooth extension of the 
normal component of v along dM to the interior of M. Subtraction of grad b(v) is 
necessary as volume preserving diffeomorphisms of a manifold with boundary leave 
the boundary invariant. 

DifFeomorphisms Leaving the Boundary Pointwise Fixed. We define the 
set 

^M,fix( M ) = i 7 ! 6 V ^( M ) I = x for all x G dM], 

that is, the volume preserving diffeomorphisms that leave the boundary pointwise 
fixed. 
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Theorem 2.1 (Ebin and Marsden (1970)) T> s ^(M) is a smooth subgroup of 
T> S (M) with Lie algebra T e D* fix (M) consisting of the space of divergence free H s 
vector fields that vanish on dM. 

In other words, this group corresponds to the choice of the no slip boundary condi- 
tions u = for the averaged Euler equations. 

We can now define P° : T V V S {M) -» r„X>" ifix (M) by 

P%(X) = [(1 - a 2 A)- 1 P e (l - a 2 A)(X o r,' 1 )] o r,. (2.1) 

Alternative Boundary Conditions. As is described in Holm, Marsden and 
Ratiu [1998a], there is an alternative choice of boundary conditions for the averaged 
Euler equations that are more like normal boundary conditions than no slip condi- 
tions. We will describe these conditions briefly and the geometry and analysis of 
this situation is explored in detail in Marsden, Ratiu and Shkoller [1999a]. 

Let n denote the outward unit normal vector field along the boundary of M and 
let N denote the corresponding normal bundle. We set 

^normalW = {r? E D'(M) | T V \ dM :N^N}. 

This group normSli \(M) corresponds to an alternative choice of boundary con- 
ditions, namely those divergence free vector fields u that are tangent to dM and 
satisfy the boundary condition 

(V n «, v) = H n {u, v) = (S n {u),v) , (2.2) 

at points of dM, for all divergence free vector fields v that are tangent to the 
boundary, where H n (or equivalently S n ) denotes the second fundamental form of 
the boundary. 



3 Mean Hydrodynamics on the Subgroup T) s ^{M) 

H l Metric on T> s fix (M ). In this section, we shall consider geodesic motion of the 
weak H 1 right invariant metric on the group fix (M) which is defined as follows. 
For X, Y € T e V^ &x (M), we set 



(X,Y} 1 = / ((X(x),Y(x)} + a 2 (VX(x),VY(x)))v(x), (3.1) 



M 



and extend (•, -)i to T> s fix (M) by right invariance. 



Submanifold Geometry. Theorem j^y ensures that T> s fix (M) is submanifold of 



V s (M) in the strong H s topology. With the induced metric ( |3.1] ), the submanifold 
P^ fix (M) also inherits weak H 1 structure. Our definition of the projection P a in 
fl2.1p from tangent spaces of V s (M) to tangent spaces of P^ fix (M) is orthogonal 
with respect to this weak H 2 structure, but P a is not a Hodge projection. In order 
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to obtain a Hodge projection which is H 1 orthogonal, one redefines the weak H 1 
metric in terms of the de-Rham Laplacian which is related to the rough Laplacian 
A = Tr VV by an additional term involving the Ricci tensor on M (see, for example, 
Section 2.3 of Shkoller [1998]). 

Using the de Rham Laplacian, Shkoller [1998] has obtained expressions for the 
unique H 1 Riemannian connection on D*(M) when M is a compact boundaryless 
manifold. Using this formula, together with submanifold geometry, he has proven 
that the weak H 1 curvature operator on P*(M) is a continuous trilinear map in 
the strong H s topology. This immediately implies the existence of unique solutions 
to the Jacobi equations along the geodesies of the H 1 Riemannian connection, and 
hence establishes the foundations for a Lagrangian stability analysis. All of this 
analysis extends trivially to manifolds with boundary on the subgroup Tt^^M) if 
one replaces the H 1 Hodge projection with our projection P a . 

Euler-Poincare equations on T e D* fix (M). The Euler-Poincare theorem (for 
background, see, for example, Marsden and Ratiu [1999] and Holm, Marsden and 
Ratiu [1998a]) uses an additional fact coming from the group structure of the prob- 
lem, namely that the right translation maps are smooth on T) 8 ^{M) which can be 
used to translate geodesic motion over the entire topological group into motion in 
the "Lie algebra" T e D* fix (M). We shall state this theorem in this context. 

Theorem 3.1 (Euler-Poincare) Equip £>^fi x (-^0 with the right invariant metric 
Then, a curve n(t) in ^^(M) is a geodesic of this metric if and only if 
u (t) = T r,(t) R r,(t)- 1 'n{t) = V(t) ° Vit)' 1 satisfies 

±u(t) = -P«o & &l {t) u(t) (3.2) 

where ad* is the formal adjoint of&d u with respect to the metric (■, -)i at the identity, 
i.e., 

(ad* v,w)i = (v, [u,w])i 

for all u,v,w eT e ^(M). 

Next, we shall prove that a unique continuously differentiable geodesic spray 
exists for the right invariant H 1 metric obtained by right translating (•, -)i over the 
entire group, but we note that even if a given metric does not have a C 1 spray from 
H s into H s , but there is still an existence theorem for geodesies, then the theorem 
is still true, when appropriately interpreted. 

We note here that use of the right invariant H 1 metric (as opposed to the usual, 
or naive, H 1 metric) is essential. For example, denote by A4 1 (A4°) the manifold of 
H s maps from T n into T n with the weak H 1 (L 2 ) topology induced by the usual H l 
(L 2 ) metric. On the orthogonal complement of Ker(l — A), M 1 is totally geodesic 
in A4°, and so the geodesic equations for both metrics are the same. 

Now, in the context of the Euler-Poincare theorem we have just stated, a straight- 
forward computation of the geodesic spray of the right invariant metric restricted 
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to right invariant vector fields on g x (M) shows that the spray coincides with the 
equations expressed in terms of the coadjoint action (of the group T) 8 qJ^M) on the 
dual of its Lie algebra) . Equivalently, there are two natural connections given by the 
Riemannian structure and right translation, and the Euler-Poincare equations can 
be obtained from subtraction of the latter from the former. For the general Euler- 
Poincare theorem for arbitrary Lagrangians, see Theorem 5.1 and the Appendix of 
Bloch, Krishnaprasad, Marsden, and Ratiu [1996] and Holm, Marsden and Ratiu 
[1998a] for the semidirect product theory. 

The form ( |3.2[) is obtained by expressing the variations 5u(t) of curves u(t) € 
T e T> s g x (M) in terms of the variations 5r](t) of curves r/(t) in T> s fix (M). In particular 
let T] : U C M? -> £>* )iix (M) be a smooth map and let 



u(t, s) 



and 



T^t^Rrj^s)- 1 



dy(t,s) 
dt 



drj(t, s) 



ds 



Then (du/ds) - (d£/dt) 



[u, £] and we obtain 
5u(t) = £ — ad u £. 



In the case of the right invariant H 1 metric on ^fi x (M), a straightforward 
computation (as in Holm, Marsden and Ratiu [1998a]) shows that 



ad* u 



a 



'A)" 1 [V tt(t)! ;(i)- a 2 [V«(f)f- Au(t)] 



where 



v = (l- a 2 A)u. 

Thus, the Euler-Poincare equation u = —P a ad* u takes the form 

,-'-^,,, M V - ,,/»> -gradp(t) 



v(t) + V u(t) v(t) - a^Vuit)} 1 ■ Au{t) 



(3.3) 



together with the divergence constraint divu = 0, the no slip boundary conditions 
it = on the boundary dM, and the initial conditions u(0) = uq. We call this 
equation the averaged Euler equations or the Euler-a equations. 

In components, and with the standard summation conventions and choice of 
Levi-Civita connection and rough Laplacian, these equations read as follows 



a 



.'/' 7 (^ + rLu*)(A»), 



-g 



ij 



dp 



where 







(A ^ = l^ + r --H^ + r 



dx l 



ml 



U 



Being Euler-Poincare equations, of course these equations share all the properties 
given by the general theory, such as a Kelvin-Noether theorem, a Lie-Poisson Hamil- 
tonian structure and so on. 
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The geodesic spray on V s fix (M). Now we are ready to state and prove the first 
main result of the paper. 

Theorem 3.2 For s > n/2 + 1, there exists a unique continuously differentiable 
geodesic spray of the metric (•, -)i on the group fix (M). 

Proof. We compute the first variation of the action function 



which we decompose as 



z JR 



and 



We have 



2 

z Jr 



£l (v) = y/ / (V(77or ? - 1 )(y),V(7?or ? - 1 )(y)) 2/ ^ 



R J M 



a 



^- I / {Vr ] {x)-[Tr 1 {x)}-\V^{x)-[Tr ] {x)]- 1 ) v{x) dxdt. 

Let e ^ r/ £ be a smooth curve in V s ^ fix (M) such that rf = r\ and (d/de)\ e= Qr] e = Sij. 
Then 



df 1 ^)-^ = a 2 



{Vrf ■ [Trf]' 1 }, Vr) • [Tr]}- 1 ) dxdt 

T)(x) 



a 





'(- 


f 

1 M 


\de 



{V^j-io^.Vi)-^]- 1 



r;(a;) 



- (V<5r/ • [Tr,]- 1 , (Vr) • [T^-^^Vt) • [Tr,]- 1 )) 



dxdt. 



Now fix x e M, and let s i— > a s be a smooth curve in M such that o"o = x, and 
(d/<is)| s= oo- s = [Tr/(x)] _1 • a for some a G T n t x \M. Define the parameterized surface 
(s,e) i ^ /(s,e) := rf oa s . Denote (d/de)(-) by (•) and (d/ds)(-) by (•)'. Then, 



and 



so that 



/(0,0) = r?(x), /(0,e) = r/ e (x), /(a, 0) = r^), 
/'(a, e) = Trf ■ a' s , f(s, e) = (d/de), rj e (a s ), 
/'(0,0)=o, and/(0,0) =5r/(x). 
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Hence, for any a £ T V ^M, and using L(a) to denote linear operation on a by the 
linear operator L, 



D 

de 



V 



D 

[T^(x)]-i(o> = J e 



_ D 


D 


o at 


o ds 



and since 



T> 


D 




e = D 


D 


de 


o ds 


r) 




ds 


o de 



7) e + i?(/,/ / )|o,or? e , 



where i? is the curvature tensor on M, we get 



^[Tr ) (x)]- 1 (a)V e — ^ 7 [Tr ! (x)]- 1 (a) 



D 

de 



r) e + R(f,f%,oV 



Thus, 



(V<5r/, (Vr/ • [Tr/]- 1 )' • (Vr) • [Tr,]- 1 ) ■ [Tr,]- 1 ' 



dxdt. 



Now, using the same argument as above, we find that 



D 



V(-5 ?? )-[T7 ? ]- 1 ,V7?-[Tr / ]- 1 



V8r] + R(Tt], r))Sr], Vr) • [Tr/] ■ [Tr]] 



= - ^V<5r/, j t {Vr) ■ [Tr]]' 1 ■ [Tr/p 1 '}^ + (i?(r), -)V(r) o r?' 1 ), <5r/> 



Finally 
D 



-^{Vr) ■ [Tr]]' 1 ■ [Tr]]' 1 '} = V-^r) ■ [Tr]]' 1 • [Tr]]' 1 ' + i?(r), Try),) ■ [Tr/]- 1 • [Tr,]- 1 ' 

-(VrHTr,]- 1 ) • (VrHTrT 1 ) • P^]- 1 ' 
-(Vr).^]- 1 ).^.^]- 1 )'.^]- 1 '. 

Integrating by parts, noting that the boundary terms vanish by virtue of the sub- 
group T>^ fix (M), and denoting by V* the L 2 formal adjoint of V, we have that 

df 1 ^)-^ = a 2 [ / -((5r/,V*{V(VM)r?- [Tr/]- 1 • [Tr/]- 1 '} 

"(VTj-tTr?]- 1 )*-^-^]- 1 )-^]- 1 ' 



+V* 



-(Vr/tTr/]- 1 ) • (Vr^fTr/]- 1 ) • [Tr/] 



i-i* 
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+(Tr R(Vfj- [Trj]" 1 , 77) • + Tr R(f), •) Vr) • [Tt/]" 1 , <5t?) 
-(V*(fl(j7,-)'7-[r77]- lt ,<yr/). 

Computing the first variation of f °, we obtain 



d£°(r,)-6r, = 11 

JR JM 



if ■ 1) ) tinll — I j ( — <$T7,») ) (/.(•(// 



/ 

^ • r \ J , 

— J), 077 ) (feat. 



Setting d£ ■ Sr/ = 0, and using the projector P a given by (|2.1|) gives 

P% o = P, a o (1 - a 2 A,)" 1 [V* {-(VviTvryiV^Tr,}- 1 ) 

+V7)[T7 ? ]- 1 V7i[T??]- 1 + (V??[Tj ? ]- 1 )(V??[T ?? ]- 1 )*} [T V ]- lt ] 
-{Ti[R(V'qTr l ~ l (-), 77) • +#(77, •) Vt?^ 1 (•)] 

+V*{ J R(7),Tr ? - 1 ') ? 7}] , (3.4) 

where 

A, = -V*[V(-)(T, ? )- 1 (T, ? )- 1 *], (3.5) 

and where again Trj^ 1 ^) defines linear operation on a vector; explicitly, operating on 
a vector V(x), this is just shorthand notation for [Tr^x)] -1 (V(x)} = [Tr](x)]~ 1 -V(x). 
Denote the quadratic form on the right-hand-side of (|3.4j ) by P^ o F°. Notice that 
F° has no derivative loss and by standard arguments is continuous from H s into 
H s . This is analogous to the equation (3.7) of Theorem 3.3 in Shkoller [1998], so we 
shall follow the proof of this theorem. The following is a combination of Lemmas 
3.1 and 3.2 in Shkoller [1998], where the detailed proofs can be found. 

Lemma 3.3 Let 

A ( .) : U,^ (M) H; (TM)|P; fix (M) U^^H^^TM^V^M) 



be defined by equation (3.5) and be the identity on P^fM). Then A(.) is a C 1 
bundle map. Also, the operator (1 — A(.)) _1 is a C 1 bundle map as well. 

We also have the following lemma which is Lemma 4 of Appendix A in Ebin and 
Marsden [1970]. 

Lemma 3.4 The L 2 Hodge projection r\ — > P^ is smooth as a function of 77. 

Choosing a local chart (U,x l ) on M, we may express the geodesic spray S : 
T^l^M) - TT ^m( M ) in thi s chart by 

Sv(v) = j t (v,v) = Mld-P^ofV^-^or?- 1 )] 
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(3.6) 



where 



T v(v, V) = <j rjfc {rf ° V X ) I '/" ° '/ ' J — f- - '/• 
The fact that the projection P° is C 1 in 77 follows from Lemma 



■ h — 1 

r] or] 



_d_ 

dx* 



and Lemma 

U since we may write Q as P%(X) = (1 - a 2 A^)- 1 P J? (l - a 2 A,,)A". 

Lemma 3.3 also shows that is C 1 as well. Hence, the term P^o (F^ — T r] (fj, fj)) 
is C 1 in T] from H s into IP. 

As for the remaining term, 

(Id - P%) o V„ ori -i(v o r?- 1 ) = (1 - a 2 A)- 1 gradA- 1 (l - c^AJdivV^-i (77 o t?" 1 ) 

+ (1 - a 2 A)- 1 gradA~ 1 [div, (1 - a 2 A)]V^ or? -i (?) o 77- 1 ). 

The fact that the first term on the right-hand-side is a C 1 map from H s into .fP 
follows identically the argument used by Ebin and Marsden [1970] to prove smooth- 
ness of the spray of the right invariant L 2 metric on D^(M). The second term is 
C 1 by similar reasoning as [div, (1 — a 2 A)] is a differential operator of order 2. This 
proves that 5 is C . 

The standard theorem for existence and uniqueness of ordinary differential equa- 
tions in a Banach space shows that unique solutions locally exist and depend 
smoothly on initial conditions. Q.E.D. 



We remark that although the statement of Lemma 3^ is C 1 , it should be clear 
that the C k result could be obtained for any positive integer k by (a messy) induc- 
tion; this type of reasoning led to the C°° regularity of the L 2 Hodge projection P 



in Lemma 3.4 



As we remarked earlier, the method of proof we have given for the smoothness of 
the spray of the right invariant H 1 metric on T> s ^ £ ix (M) also provides well-posedness 
for the one dimensional Camassa-Holm equation on S , since the spray in this case 
does not have a loss of derivative (c.f. Remark 3.5 in Shkoller [1998]). The situation 
of the present paper is, however, complicated in a nontrival way by the presence of 
boundaries. We state the situation for the CH equation as the following theorem. 

Theorem 3.5 The Cauchy problem for the ID CH equation J\ ), given by 



V 



with initial conditions 



+ - 2 (v 



T] 



-1)2 
hi 



V 



(3.7) 



V (0) 



77(0) = n , 



has a unique solution (77,77) in V 8 ^ 1 ) x H 8 ^ 1 ) for s > § on a finite time interval 
where the solution has C 1 dependence on time and smooth dependence on initial 
data. 



Theorem 3.2 implies many other results too, such as finite smoothness with re- 
spect to initial conditions in Eulerian representation if one is willing to map between 
different Sobolev spaces, smoothness with respect to the time variable, etc. 
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Remark. The subgroup fix (M) ensures that the boundary terms vanish in the 
above variational principle, as the variation 5rj(x) = for all x £ dM. For variations 
which are not required to vanish on the boundary, one obtains natural boundary 
conditions on dM, which are discussed in Marsden, Ratiu and Shkoller [1999a]. 



4 The limit of zero viscosity on V s ^(M) 

The Navier-Stokes a-model is obtained by adding viscous diffusion to the Euler-a 
model. The equations are given by 

d t u - vAu + (1 - a 2 A)' 1 [V u (l - a 2 A)u - a 2 Vu* ■ Au] = -(1 - a 2 A)' 1 grad p. 

(4.1) 

Foias, Holm, and Titi [1999] have proven global well-posedness of ( |1.1| ) and have 
obtained estimates on the dimension of the global attractor. Chen et al. [1999] have 
given some very interesting results concerning the applicability of these equations to 
high Reynolds number flows in channels, including relations to the law of the wall. 

Having proven the smoothness of the geodesic spray of the Euler-a equations, 
we follow Ebin and Marsden [1970] and use a product formula approach to prove the 
existence of the limit of zero viscosity of (|4.1| ). In the case that a = 0, this limiting 
procedure is believed to be valid only for compact manifolds without boundary 
(e.g., for flows with periodic boundary conditions), as the Navier-Stokes equations 
and the Euler equations do not share the same boundary conditions on manifolds 
with boundary. When, a / 0, however, as we shall discuss in the last section, a 
certain type of elasticity is added into the Euler-a model, and the mean motion 
of the fluid exhibits normal stress effects. Because of this, we may prescribe zero 
velocity boundary conditions even in the inviscid limit, and thus extend the limit of 
zero viscosity theorems for the averaged Euler equations to manifolds with boundary. 

The following is the Euler-a version of Theorem 13.1 of Ebin and Marsden [1970]. 



Theorem 4.1 Let S : TV S -(M) -> TTV S „{M) be the Euler-a vector field given 
by Q). Let T : T e ^ fix (M) -» T e t> s - a {M) be a given map, where a > 2. Let 
T be a bounded linear map which generates a strongly- continuous semi-group Ft : 
T e V 8 zjM) -> T e D*(M), t > 0, and satisfies \\F t \\ s < e 131 for some (3 > and some 
s. Extend F t to TP^ fix (M) by 

Ft\X v ) = TR n ■ F t ■ TR r) -i(X rj ) 

for X v G T V V S ^ M {M), and let f be the vector field f : TV S ^{M) -> TTV s - a {M) 

associated to the flow F t . 

Then S + vT generates a unique local uniformly Lipschitz flow on TD* fix (M) 
for v > 0, and the integral curves c u {t) with c"(0) = x extend for a fixed time r > 
independent of v and are unique. Further, 

lim c»(t) = c°(t) 

u— »0 

for each t, < t < r, the limit being in the H s topology, s > (n/2) + 1 + 2a. 
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See Lemmas 2-5 in Appendix B of Ebin and Marsden [1970] for an explanation 
of 2a in the hypothesis of Theorem 4.1. Since Theorem |3.2| proves that the geodesic 



spray is smooth on "D 8 ^{M), we must verify the hypothesis of Theorem 4.1 for 



-A. 



Lemma 4.2 The operator A generates a strongly continuous contraction semi- 
group on the space T e V s ^ fix (M), for s > (n/2) + 1. 

Proof. Standard energy estimates give the result. See, for example, Theorem 3.2 
of Temam [1988], page 70, and Theorem 13.2 of Ebin and Marsden [1970]. Q.E.D. 

Remark. With initial data uq € T e X>~ fix (M), the solution u(t) is also C°° as 
consequence of the regularization of parabolic flows. 

The way one proves the product formulas mentioned in the introduction is ex- 
actly the same as is outlined in Ebin and Marsden [1970] and Marsden, Ebin and 
Fischer [1972]; one uses the general theorems on product formulas (formulas of Lie- 
Trotter type) for smooth vector fields on infinite dimensional manifolds, as discussed 
in, for example, Chorin, Hughes, Marsden and McCracken [1978] and in Abraham, 
Marsden and Ratiu [1988] using the smoothness of the spray and the projection 
operator. For product formulas involving the Stokes operator, one couples the 
smoothness of the spray with the analytic property of the generated semigroup, 
as explained in Ebin and Marsden [1970]. Of course this technique is closely related 



to Theorem 4.1 discussed above. 



5 Second-grade Fluids 

As noted in Chen et al. [1999], non-Newtonian fluids of second grade satisfy the 
same type of equations as the mean motion model equations for incompressible 
ideal fluids. In fact, we shall show that in the incompressible case, the equations are 
identical in the limit of zero viscosity. For simplicity we shall make the presentation 
in Euclidean space. (See Marsden and Hughes [1983] for how to put this discussion 
on manifolds). 

The constitutive equation of fluids for a fluid of grade two is given by (see, for 
example, Noll and Truesdell [1965]) 

T = -pl + vA l + a l A 2 + a 2 Al, (5.1) 

where T is the stress tensor, p is the hydrostatic pressure, 1 is a unit tensor, v is 
viscosity, a.\ and a 2 are normal stress moduli. Here, Ai and A2 are Rivlin-Ericksen 
tensors (Rivlin and Ericksen [1955]) defined by 

Ai = Vu + (Vu)* , 
DAi 

A 2 = — - + (V«)*Ai + AiVn . 
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Here, the tensors are of type (1,1) and D/dt = (d/dt) + u • V is the material 
derivative. 

To satisfy the Clausius-Duhem inequality, the first material constant has to 
satisfy v > and the constants a\ and «2 satisfy a\ + a 2 = . If in addition 
one asks that the specific Helmholtz free energy is a minimum when the fluid is in 
equilibrium, then 

ai > 0. 

The results expressed by the last inequalities are derived in Dunn and Fosdick [1974]. 

We now turn to the incompressible case, i.e., divu = 0. Under this assumption, 
we make the following calculations on the stress tensor: divT : 

divT = - div pi + v div Ai + «i(div A 2 — div Af) , 

divAi = div(Vu)* = Au , 
d 

divA 2 = — div Ai + div((u • V)Ai) 

+ div(AiVu) + div((Vu)* A 1 ) , 
div((n-V)Ai) = Vu ■ VA X + u ■ V(Au) , 

div(AiVu) = J](Au) fe Vu fc + ^V(Tr(Vu) 2 ) + ^V(TY(Vu(Vu)')) , 

k 

div((Vn)*Ai) = Au • Ai + (Vu) ( • VAi , 

div(A?) = An • Ai + (Vu)* • VAi + Vu ■ VA 1 , 
d 

div A 2 - div(A?) = -(Au) + u ■ V(Au) 

+(Vu)* • An + ^V(Tr(Vn) 2 ) + -V(Tr(Vu(Vu)')) . 



Substituting div T into the balance of momentum with no body force acting on the 
fluid, 

Du ,. m 

"iir = dlTT ' 

and making use of the calculations above, we obtain the following equation for a 
second-grade fluid 



■(1 - «i A)u = i)Au + ai (Vu)' • Au + V «i D : D - - , (5.2) 
dt V P, 



where D: D = Tr(DD*) and 



OL\ = , V = - 

P P 



and the stretch tensor is given by 



D = -Ai. 

2 



18 



The theory of second-grade fluids has been developed over the last 40 years by 
a number of people, such as Noll and Truesdell [1965], Coleman and Markovitz 
[1964], Markovitz and Coleman [1964], and Coleman, Markovitz, and Noll [1966]. 
These authors have shown that second grade fluids can be regarded (like the Navier- 
Stokes theory) as an approximation to the general theory of simple fluids with fading 
memory, valid in sufficiently slow flows. 



Inviscid Second Grade Fluids. Letting the viscosity v = in (5.2) and defining 
o? = a\ , v = u — a 2 An gives the following equation 

Dv 

— a 2 (Vu)' • An = — grad P , divu = divv = , (5-3) 

where the second-grade fluid pressure function is given by 

-2 2 

P = 1 _ 2L Tr(Vn • (Vu)*) - — TY(Vu) 2 , 
p 2 2 

and where P is the generalized pressure. 



We observe that the equations (5.c) are identical to the averaged Euler equations 
(3.5). In either case, one can eliminate the pressure terms in the standard fashion 
using the projection to the divergence free part or, if one prefers, one can view the 
pressure as being determined implicitly from the incompressibility constraint. While 
these equations are the same, the physical interpretation and the basic mechanics 
are different. 



Comments on the Limit of Zero Viscosity. The equations for a second grade 
fluid with viscosity ( |5.2| ) are not the same as the Navier-Stokes— a equations ( |l.3|) . 
The difference is simply that in the case of second grade fluids, one uses the Lapla- 
cian of the velocity vector u while in the Navier-Stokes— a equations one uses the 
Laplacian of the momentum v. 

The theory for second grade fluids is, in a sense, much easier than that for the 
averaged Navier-Stokes equations. This is because when the equations are written 
with the dynamical variable u, there is no derivative loss in the dissipative term and 
so the vector field in Lagrangian coordinates is already smooth. The result on the 
limit of zero viscosity for second grade fluids corresponding to that for the Navier- 
Stokes— a equations is elementary, following simply from the smoothness of the flow 
of a smooth vector field as a function of parameters. 

In the literature on second-grade fluids it is generally assumed that the viscos- 
ity does not vanish. However following the approximation procedure of Coleman 
and Markovitz [1964], and Noll and Truesdell [1965] one can show that the case of 
zero viscosity does not contradict their theory. Hence for a simple fluid with fading 
memory, a second-grade fluid with zero viscosity furnishes a complete second-order 
correction, for normal stress effects, to the theory of perfect (in the terminology of 
Noll and Truesdell) elastic fluids, in the limit of slow motion. Notice that ci\ is a ma- 
terial constant in second-grade fluid theory, whereas a is the fluctuation amplitude 
scale parameter (i.e., a flow regime parameter) in the average Euler equations. 
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Another aspect to be mentioned here is that the theory of simple fluids obeys 
the basic mechanical principles: the principle of material frame-indifference, which 
assert that the response of a material is the same for all observers, and the principles 
of determinism and local action, which asserts that the present stress at a particle is 
determined by the history of an arbitrary small neighborhood of that particle. The 
same principles are true of the averaged Euler equations and the Navier-Stokes— a 
equations. 



Existence and Uniqueness. In the viscous case, existence and uniqueness has 
already been shown by Cioranescu and Gorailt [1997] and Galdi [1995], although as 
we have mentioned, many results, such as short time existence and uniqueness and 
the limit of zero viscosity results follow from our work too. 

The question of the existence and dimension of attractors for the Navier-Stokes— a 
equations is addressed in Foias, Holm and Titi [1999]. We note that these techniques 
do not allow one to take the limit of zero viscosity, whereas ours do allow such a 
limit. 



Specific Solutions. Several formal solutions (ignoring boundary conditions) have 
been given for planar homogeneous incompressible flows of second grade fluids. The 
stream function formulation of the dynamics of such flows in two dimensions, with 
stream function ip and velocity u related by 

\ dy ' dx ' / ' 

is given by 

qt + V„g = i/Aw, (5.4) 

where 

g = -(l-a 2 A)A^ and uj = -Atp. 
iFor example, Ting [1963] provided a set of exact solutions for planar startup flows. 



Further, Rajagopal and Gupta [1981] showed that equation (5.4) has separable so- 
lutions that may be interpreted damped lattice of vortices, two impinging 
rotational flows, and a viscoelastic analog to Kelvin's "Cat's eyes" vortices. These 
formal separable solutions impose the condition 

Ai/j = cip, (5.5) 

where c is a nonzero constant. Of course, this condition implies the linear relation 
q = —(1 — a 2 c)cxp. Hence, the transport term V u q = J(q,ip) vanishes in flO) ) and 
the remaining terms impose 

vq 

Qt — z g - , 

1 — ore 
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which implies exponential time dependence for q and, thus, for tp. For example, the 
stream function for the viscoelastic analog to Kelvin's "Cat's eye" vortices found in 
Rajagopal and Gupta [1981] is given in separated form by 

ip = A cosh(ax) cos(6y) exp(Ai) . (5-6) 



This stream function is a solution of (|5.4| ) under condition ( |5.5|) for arbitrary real 
A,a,b, with c = a 2 — b 2 and A = vj(\ — a 2 c). Relation (|5.6|) with A = and 



A,a,b arbitrary is also a stationary solution of equation (5.4) in the limiting case 
that v = 0. 

Two-dimensional incompressible Navier-Stokes-a dynamics may be expressed in 
a form similar to equation ( |5.4|) , namely, 

qt + y u q = vAq, (5.7) 

where 

q = z ■ curly = —(1 — a 2 A)Ai/> 
is potential vorticity. Thus, under condition (|5.5|) (and still ignoring boundary con- 



ditions), separable solutions for second-grade fluids satisfying (5.4) become separable 
solutions for the two-dimensional Navier-Stokes-a model satisfying ( |5.7D under the 
change u — *■ (1 — a 2 c)z/. In the case v = 0, formal solutions satisfying the lin- 
ear relation ( |5.5| ) are common to both theories. Another example of a steady flow 
that is common to both theories in the case v = is the Gaussian jet, for which 
ip = exp(-y 2 ). 

For the ideal limiting case in which v = it should be clear that steady incom- 
pressible flows of both second grade fluids and the two-dimensional Euler-a model 
occur whenever q and ip are functionally (not just linearly) related. For example, 
for periodic boundary conditions, critical points of the functional 



H* = H + C$ = i y [# + $(q)]dxdy , 



(the sum of the energy and the domain integral of an arbitrary function $ of the 
potential vorticity) are steady flows of both models. This observation leads imme- 
diately to a unified viewpoint of nonlinear stability criteria for steady periodic flows 
of either model, following the energy-Casimir method, as in Holm, Marsden, Ratiu 
and Weinstein [1985]. 

Exact solutions of the three-dimensional Navier-Stokes-a model may also be 
sought by adapting the classical solutions for the two-dimensional Navier-Stokes 
equations, e.g., for flows between two rotating disks, Beltrami flows, etc., as re- 
viewed, e.g., in Wang [1991]. 
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